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We generalize the total angular momentum method for computing Cosmic Microwave Background 
anisotropies to Friedman-Robertson- Walker (FRW) spaces with arbitrary geometries. This unifies 
the treatment of temperature and polarization anisotropies generated by scalar, vector and tensor 
perturbations of the fluid, seed, or a scalar field, in a universe with constant comoving curvature. 
The resulting formalism generalizes and simplifies the calculation of anisotropies and, in its integral 

ON ■ form, allows for a fast calculation of model predictions in linear theory for any FRW metric. 
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>y I. INTRODUCTION 

a ■ 

' The study of the Cosmic Microwave Background (CMB) radiation holds the key to understanding the seeds of the 
\q . structure we see around us in the universe, and could potentially enable precision measures for most of the important 
fSJ ' cosmological parameters. For this reason, as well as because of its intrinsic interest, one would like a physically 
transparent framework for the study of CMB anisotropies which is as general, powerful, and flexible as possible. 

Theoretically, the calculation of CMB anisotropies is "clean" , involving as it does only linear perturbation theory. 
However the calculations can become quite complex once one allows for the possibility of non-flat universes, non-scalar 
perturbations to the metric, and polarization as well as temperature anisotropies. Recently Hu & White |l[ presented 
a formalism for calculating CMB anisotropies which treats all types of perturbations, temperature and polarization 
ON ' anisotropies, and hierarchy and integral solutions on an equal footing. The formalism, named the total angular 
momentum method, greatly simplifies the physical interpretation of the equations and the form of their solutions (see 
e.g. [||). However it was presented in detail only for the case of flat spatial hypersurfaces. Here we generalize the 
treatment for the curved spaces of open and closed Friedman-Robertson- Walker (FRW) universes. 
' Aspects of this method in open (hyperbolic, negatively curved) geometries have been introduced in Hu & White 
^p" 1 H and Zaldarriaga, Seljak & Bertschinger [Q for the cases of tensor temperature and scalar polarization respectively. 
The latter work also addressed methods for efficient implementation through the line of sight integration technique 
In this paper, we complete the total angular momentum method for arbitrary perturbation type and FRW metric, 
paying particular attention to the case of open universes because of its strong observational motivation. As an example 
we use this formalism to compute the temperature and polarization angular power spectra of both scalar and tensor 
modes in critical density and open inflationary models. We incorporated the formalism into the CMBFAST code of 
, Seljak & Zaldarriaga ||, which has been made publically available. 

The outline of the paper is as follows: we begin by establishing our notation for fluctuations about a FRW back- 
ground cosmology in §|n| We then present the Boltzmann equation in our formalism in §[H, which contains the main 
results. We give some examples and discuss applications in §|y|. Some of the more technical parts of the derivations 
(the Einstein, radial and hierarchy equations) are presented in a series of three Appendices. 
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II. METRIC AND STRESS-ENERGY PERTURBATIONS 



In this section, we discuss the representation of the perturbations for the cosmological fluids and the geometry of 
space-time. We start by defining the basis in which we shall expand such perturbations and their representation under 
various gauge choices. 

We assume that the background is described by an FRW metric — a 2 7 Ml/ with scale factor a(t) and constant 
comoving curvature K = —Hq(1 — iltot) in the spatial metric 7^. Here greek indices run from to 3 while latin 
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indices run over the spatial part of the metric: i, j = 1, 2, 3. It is often convenient to represent the metric in spherical 
coordinates where 

lid dx l dx3 = IK]- 1 [d X 2 + sm 2 K X (d6 2 + sin 2 d(j> 2 )] , (1) 

with 

sin* ( X ) = { Sin ?W ' £<°' (2) 

where the flat-limit expressions are regained as K — > from above or below. The component corresponding to 
conformal time 

is 7oo = -1. 

Small perturbations h^ v around this FRW metric 

9fiu = a 2 {~f^ + hftv) , (4) 

can be decomposed into scalar (m = 0, compressional) , vector (m = ±1, vortical) and tensor (m = ±2, gravitational 
wave) components from their transformation properties under spatial rotations 



A. Eigenmodes 

In linear theory, each eigenmode of the Laplacian for the perturbation evolves independently, and so it is useful to 
decompose the perturbations via the eigentensor Q^ m \ where 

V 2 Q('«) = - -k 2 Q im \ (5) 

with "|" representing covariant differentiation with respect to the three metric jij. Note that the eigentensor Q( m ) 
has \m\ indices (suppressed in the above). Vector and tensor modes also satisfy the auxiliary conditions 



)(±i)|i 



0. 



7 «Qg*) =Q (±3)|i = 0j ( 6 ) 

which represent the divergenceless and transverse-traceless conditions respectively, as appropriate for vorticity and 
gravity waves. In flat space, these modes are particularly simple and may be expressed as 

Qit™i m (^i ^ ^2)11 • • • (ei ± «e 2 )i m exp(ifc • x) , (K = 0, m > 0) , (7) 

where the presence of ej, which forms a local orthonormal basis with £3 = k, ensures the divergenceless and transverse- 
traceless conditions. 

It is also useful to construct (auxiliary) vector and tensor objects out of the fundamental scalar and vector modes 
through covariant differentiation 

Q? ] = -k-'Qf , Qf = k- 2 Q^ + ^Q (0) , (8) 

n (±l) _ _/„n-l/«(±l), n (±l) 



- 'J 



(2*)- 1 (Qg 1J +Q}J 1 0. (9) 



The completeness properties of these eigenmodes are discussed in detail in H, where it is shown that in terms of 
the generalized wavenumber 



q=y/k* + (\m\ + l)K, v = q/\K\, (10) 

the spectrum is complete for 



v > 0, K < 0, 

= 3,4,5..., K>{). 



(11) 



A deceptive aspect of this labelling is that for an open universe the characteristic scale of the structure in a mode 
is 27r//c and not Zit/q, so all functions have structure only out to the curvature scale even as q — » 0. We often go 
between the variable sets (k,rf), (q,rj) and (v,x) f° r convenience. 
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B. Perturbation Representation 



A general metric perturbation can be broken up into the normal modes of scalar (to = 0), vector (to = ±1) and 
tensor (to = ±2) types, 

m 
m 

h l3 = £ 2 J ffi m) Q (m) 7y + . (12) 

m 

Note that scalar quantities cannot be formed from vector and tensor modes so that = and H^ 1 ' — for to ^ 0; 
likewise vector quantities cannot be formed from tensor modes so that £?( m ) = for |to| = 2. 

There remains gauge freedom associated with the coordinate choice for the metric perturbations (see Appendix 
A 2 ) . It is typically employed to eliminate two out of four of these quantities for scalar perturbations and one of the 
two for vector perturbations. The metric is thus specified by four quantities. Two popular choices are the synchronous 
gauge, where 

= Hl, Hr^ = hr, 
H^, ^ = hy , H^p ^ = H, (13) 

and the generalized (or conformal) Newtonian gauge, where 

A< > = *, flW=V, 

Hf } = $ , =H. (14) 

Here and below, when only the to > expressions are displayed, the to < expressions should be taken to be identical 
unless otherwise specified. 

The stress energy tensor can likewise be broken up into scalar, vector, and tensor contributions. Furthermore 
one can separate fluid (/) contributions and seed (s) contributions. The latter is distinguished by the fact that the 
net effect can be viewed as a perturbation to the background. Specifically T^ v — T^ v + 5T^ V where T° = —pf, 
T° = T * = and = ■p s b i ] is given by the fluid alone. The fluctuations can be decomposed into the normal modes 
ofSflAlas 



<*n = -E m [pA m) +Ps]Q {m \ 

= -£J(p/+p/)4 m) +^]Q (m)i , 
ST) = E m [^/ m) +PsW j Q {m) + b/4 m) +Ps)Q (m)i j - 

Since Sl* 1 ^ — Sp™^ = for to ^ 0, we hereafter drop the superscript from these quantities. 
A minimally coupled scalar field tp with Lagrangian 



(15) 
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can be treated in the same way with the associations 
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£ = —V^W^d^V + 2F(^)] (16) 



2V = |a- 2 2 + V, (17) 

for the background density and pressure. The fluctuations ip = <\> + 5(f) are related to the fluid quantities as R] 

SP4, = 5p<k + 2V,*ty = a- 2 (<j)8q) - A^tf) + V ^ , 
(P0+^)(^ O) -S (O) ) = a- 2 ^, 

P**i o) =0. (18) 
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The evolution of the matter and metric perturbations follows from the Einstein equations G M „ — SnGT^ and 
encorporates the continuity and Euler equations through the implied energy-momentum conservation T^".^ = 0. We 
give these relations explicitly for the scalar, vector and tensor perturbations in both Newtonian and synchronous 
gauge in Appendix [A| (see also ||). 

These equations hold equally well for relativistic matter such as the CMB photons and the neutrinos. However 
in that case they do not represent a closed system of equations (the equation of motion of the anisotropic stress 
perturbations is unspecified) and do not account for the higher moments of the distribution or for momentum 
exchange between different particle species. To include these effects, we require the Boltzmann equation which 
describes the evolution of the full distribution function under collisional processes. 



III. BOLTZMANN EQUATION 

The Boltzmann equation describes the evolution in time (ry) of the spatial (a?) and angular (n) distribution of the 
radiation under gravity and scattering processes. In the notation of O, it can be written implicitly as 

^-f( v , x, n) = ^T + n% = C[f] + G[h^} , (19) 

where T — (9, Q + ill, Q — ill) encapsulates the perturbation to the temperature 9 = AT/T and the polarization 
(Stokes Q and U parameters) in units of the temperature fluctuation. The term C accounts for collisions, here 
Compton scattering of the photons with the electrons, while the term G accounts for gravitational redshifts. 



A. Metric and Scattering Sources 

The gravitational term G is easily evaluated from the Euler-Lagrange equations for the motion of a massless particle 
in the background given by |p|,p|,^0[ : 



G[h^] = -rtnPhij + n l h 0i + -n*/i o|i ,0,0 



(20) 



Note that gravitational redshifts affect different polarization states alike. As should be expected, the modification 
from the flat space case involves the replacement of ordinary spatial derivatives with covariant ones. 

The Compton scattering term C was derived in in the total angular momentum language. Though the basic 
result has long been known [ [ll"|]l1| , this representation has the virtue of explicitly showing that complications due to 
the angular and polarization dependence of Compton scattering come simply through the quadrupole moments of the 
distribution. Here 



C[T] 



T(n) 



dn' 

47T 



9' + h ■ v B , 0, 



+ .L / dn' J~] P (m) (n,n')!T(ri'), 



2 



m=-2 



(21) 



where the differential cross section for Compton scattering is t = n e crra where n e is the free electron number density 
and ax is the Thomson cross section. The bracketed term in the collision integral describes the isotropization of the 
photons in the rest frame of the electrons. The last term accounts for the angular and polarization dependence of the 
scattering with 



p(m) 



V 



ml \frr 
2 I 1 2 



3 \sml \rm 
22*2 1 2 



'2 ± 2 2 1 2 



3 2 Y 2 m '. 



-2 X 2 



3 \rm/ \^ m 
2-2*2 I 2 



'-2 1 2 2 A 2 



°-2 I 2 -2*2 



(22) 



where Y™' = Y l m *(h') and s Y( al = s Y7™*(n') and the unprimed harmonics have argument n. Here s 5^ m are the 
spin- weighted spherical harmonics [p^| -pM . 
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B. Normal Modes 



The temperature and polarization distributions are functions of the position x and the direction of propagation 
of the photons n. They can be expanded in modes which account for both the local angular and spatial variations: 
s G?(x,n),i.e. 



<d(r],x,n) 
(Q±iU)(r,,x,n) 



J v ' l m=-2 

/tS^E E(^ (m) ±^ (m) ) ±2 Gr, 

J (Z7T) l m— ■> 



(23) 



with spin s = describing the temperature fluctuation and s = ±2 describing the polarization tensor. Ei and -B; 
are the angular moments of the electric and magnetic polarization components. It is apparent that the effects of the 
local scattering process C is most simply evaluated in a representation where the separation of the local angular and 
spatial distribution is explicit I3J|, with the former being an expansion in s 3j m - The subtlety lies in relating the local 
basis at two different coordinate points, say the last scattering event and the observer. 

In flat space, the representation is straightforward since the parallel transport of the angular basis in space is trivial. 
The result is a product of spin-weighted harmonics for the local angular dependence and plane waves for the spatial 
dependence: 



s GT(x,n) = H) i 



47T 

21 + V 



s y ; m (n)]exp(ifc-f) 



(K = 0) 



(24) 



Here we seek a similar construction in an curved geometry. We will see that this construction greatly simplifies 
the scalar harmonic treatment of p6| , p^7| j4[ and extends it to vector and tensor temperature || modes as well as all 
polarization modes. 

To generalize these modes to the curved geometry, we wish to repl ace the plane wave with some spatially dependent 
phase factor exp[i6(x, k)] related to the eigenfunctions Q( m ) of § 11 A while keeping the same local angular dependence 
(see Eq. C2). By virtue of this requirement, the Compton scattering terms, which involve only the local angular 
dependence, retain the same form as in flat space. In Appendix [cj we derive S G™ by recursion from covariant 
contractions of the fundamental basis Q( m >. The result is a recursive definition of the basis 



71 K s^l )\i — 21+1 L s 1 U^i-lJ _ s h l + l\s U l+l!\ ~ + j\ s U l 



constructed from the lowest /-mode of Eq. (B2) with the coupling coefficient 



(I 2 - m 2 )(P - s 2 ) 



I 2 



I 2 

I 7tK 

q 2 



(25) 



(26) 



The structure of this relation is readily apparent. The recursion relation expresses the addition of angular momentum 
and is the defining equation in the total angular momentum method. It says the "total" local angular dependence 
at (say) the origin is the sum of the local angular dependence at distant points ( "spin" angular momentum) plus the 
angular variations induced by the spatial dependence of the mode ( "orbital" angular momentum) . 

The recursion relation represents the addition of angular momentum for the case of an infinitesimal spatial sepa- 
ration. Here the leading order spatial variation is the gradient [n l ( s G™)u] term which has an angular structure of a 



Y]° with the intrinsic 



dipolc Yi . The first term on the rhs of equation (12q) arises from the Clebsch-Gordan relation that couples the orbital 



YJ m to form I ± I states, 



v/(2Z + l)(2Z-l) 



y/(2l + l)(2Z + 3) 



1(1 + 1) 



(27) 



where the coupling coefficient is s c™ = yj (I 2 — m 2 )(l 2 — s 2 )/l 



The second term on the rhs of the coupling equation (26) accounts for geodesic deviation factors in the conversion of 
spatial structure into orbital angular momentum. Consider first a closed universe with radius of curvature 1Z = K^ 1 / 2 . 
Suppressing one spatial coordinate, we can analyze the problem as geometry on the 2-sphere with the observer situated 
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at the pole. Light travels on radial geodesies or great circles of fixed longitude. A physical scale A at fixed latitude 
(given by the polar angle \) subtends an angle a = X/TZ sin \. I n the small angle approximation, a Euclidean analysis 
would infer a distance related by 

V{d) = "Ksinx = K~ 1/2 sm X , (K > 0), (28) 
called here the angular diameter distance. For negatively curved or open universes, a similar analysis implies 

V{d) = \K\- 1/2 sinhx , (K < 0). (29) 
Thus the angular scale corresponding to an eigenmode of wavelength A is 

A 1 



1Z sinh x v sinh \ 

For an infinitesimal change x, orbital angular momentum of order I is stimulated when 

r,*~[l + 0(l 2 K/q 2 )], 



(30) 



(31) 



which explains the factors of l 2 K/q 2 in the coupling term in a curved geometry. We shall see in §111 D that these 
infinitesimal additions of angular momentum and geodesic deviation may be encorporated into a single step by finding 
the integral solutions to the coupling equation (Eq). 



C. Evolution Equations 



It is now straightforward to rewrite the Boltzmann equation ( p~S| ) as the evolution equations for the amplitudes of 
the normal modes of the temperature and polarization = (0| ; , E^"" 1 , B\ f ) . The gravitational sources and 

scattering sources of these equations follow from Eq. (|2^) and ( pl| ) by noting that the spin harmonics are orthogonal, 



d£l{ a Yr){.YF*) = 6w6 v 



(32) 



The term n l Tji is evaluated by use of the coupling relation Eq. ( p5| ) for n l { s G™)|i. It represents the fact that spatial 
gradients in the distribution become orbital angular momentum as the radiation streams along its trajectory x(n). For 
example, a temperature variation on a distant surface surrounding the observer appears as an anisotropy on the sky. 
This process then simply reflects a projection relation that relates distant sources to present day local anisotropies. 
With these considerations, the temperature fluctuation evolves as 
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and the polarization as 



(m) 



(21 - 1) 



o( m ) _ 



(21 



1-1 fi (m) 
3) U;+1 



^(™) 



(l > m), 



E, 



j,(m) _ 



2^1 



(21 - 1) 



&l-i ~ 



2m 



1(1 + 1) 



D 



(m) 



B l - <I 



2^1 



-B 



(m) 



2m 



-E 



(m) 



2^1+1 F (m) 

(21 + 3) 1+1 
-B\ 



f[E[ m) + V6P (m) S L2 ] 



(21-1) 1(1 + 1) 1 (21 + 3) 1+1 

The temperature fluctuation sources in Newtonian gauge are 



-+B 



(m) 



(33) 



(34) 



sP=tvU + V, 



4 0) 



s. 



(1) 



(35) 



5. 



(2) 



tPW - H . 
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and in synchronous gauge, 



(0) 
B ' 


4 0) 


= fp(°) 


(1) 
B ' 




= fP« 






= fpW 



S{ x) = f«W , 4 1} = f P« - ^ V 1 - 2K /k 2 h v , (36) 



The I = m = 2 source doesn't contain a curvature factor because we have recursively defined the basis functions in 
terms of the lowest member, which is I = 2 in this case. In the above 

pM = _L [ei m) - \/6£4 m) ] (37) 
and note that the photon density and velocities are related to the I = 0, 1 moments as 

8 1 = 46< 0) , 4 m) = e[ m) ; (38) 
whereas the anisotropic stresses are given by 

irMggO =l 2 y^ Kn, - i^)0W, (39) 

which relates them to the quadrupole moments (Z = 2) as 

(1 - ZK/k*) 1 /**!? = , (1 - 2X/fc 2 ) 1 /^(D = Trf = ^6< 2) . (40) 

O DO 

The evolution of the metric and matter sources are given in Appendices |A 3| — |A | . 

D. Integral Solutions 

The Boltzmann equations have formal integral solutions that are simple to write down. The hierarchy equations 
for the temperature distribution Eq. ( |33| ) merely express the projection of the various plane wave temperature sources 

s i" l) o G T on the sk y toda y ( see E q- &)■ Likewise Eq. @ expresses the projection of -V6P (m) f er T ±2 G\ n . 

The projection is obtained by extracting the total angular dependence of the mode from its decomposition in 
spherical coordinates: i.e. into radial functions times spin harmonics S Y" 1 . We discuss their explicit construction in 
Appendix |b|. The full solution immediately follows by integrating the projected source over the radial coordinate, 







21 


+ 1 


E (m) 


(vo,q) 


21 


+ 1 


D l 




21 


+ i 



rVo 

/ d V e-^S : 
Jo 



(m) Ajm) 

3 





r-VO 



d V Te- T (-V6P {m) )[3l m \ (41) 

where the arguments of the radial functions (0;,e;,/3;) are the distance to the source x = V—K(rjo — if) and the 
reduced wavenumber v = q/\/—K (see Appendix^ for explicit forms). 

The interpretation of these equations is also readily apparent from their form and construction. The decomposition 
of S G™ into radial and spherical parts encapsulates the summation of spin and orbital angular momentum as well 



as the geodesic deviation factors described in §IIIB. The difference between the integral solution and the differential 
form is that in the former case the coupling is performed in one step from the source at time r\ and distance xiv) to 
the present, while in the latter the power is steadily transferred to higher I as the time advances. 

Take the flat space case. The intrinsic local angular momentum at the point (x, n) is s Yp but must be added to the 
orbital angular momentum from the plane wave which can be expanded in terms of jiY, . The result is a sum of \l — j\ 
to I + j angular momentum states with weights given by Clebsch-Gordan coefficients. Alternately a state of definite 
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angular momentum involves a sum over the same range in the spherical Bessel function. These linear combinations 
of Bessel functions are exactly the radial functions in Eq. ( fil] ) for the flat limit H . 

For an open geometry, the same analysis follows save that the spherical Bessel function must be replaced by a 
hyperspherical Bessel function (also called ultra-spherical Bessel functions) in the manner described in Appen dix [B . 



The qualitative aspect of this modification is clear from considering the angular diameter distance arguments of § III B . 
The peak in the Bessel function picks out the angle which a scale k~ l s» \/— Kv" 1 subtends at distance d « x/v~~K- 
A spherical Bessel function peaks when its argument hd w I or X/d w 9 in the small angle approximation. The 
hyperspherical Bessel function peaks at kT> = ^sinhx ~ I for v > 1 or A/2? w 9 in the small angle approximation. 
The main effect of spatial curvature is simply to shift features in Z-space with the angular diameter distance, i.e. to 
higher / or smaller angles in open universes. Similar arguments hold for closed geometries [ fL7| . By virtue of this fact 
the division of polarization into E and -B-modes remains the same as that in flat space. More specifically, for a single 
mode the ratio in power is given by 

EM m) ? r °' m = ' 



6, m = ±1, (42) 



E,[fcH a 1 8/13, m = ±2 

at fixed source distance y with z/ skiff x 3> 1. 

The integral solutions (|4l| ) are the basis of the "line of sight" method [^|Jl^1 for rapid numerical calculation of CMB 
spectra, which has been employed in CMB FAST. The numerical implementation of equations ( [il|) requires an efficient 
way of calculating the radial functions (<fii, e;, 0i). This is best done acting the derivatives of the hyperspherical Bessel 



function in the radial equations (B3)-(Bq) and (Bll) on the sources through integration by parts. The remaining 
integrals can be efficiently calculated with the techniques of Q for generating hyperspherical Bessel functions. The 
tensor CMBFAST code has now been modified to use the formalism described in this paper and the results have been 



cross-checked against solutions of the Boltzmann hierarchy equations (|33j)-(|3j) with very good agreement. 

E. Power Spectra 

The final step in calculating the anisotropy spectra is to integrate over the A-modes. The power spectra of tem- 
perature and polarization anisotropies today are defined as, e.g. Cf e = (|a im | 2 ) for 9 = £ a lm Y{ n with the average 
being over the (21 + 1) m-values. In terms of the moments of the previous section, 

2 

(21 + ifC?* = - I ^ Y q 3 X[ m) *x[ m) , (43) 

where X takes on the values Q, E and B for the temperature, electric polarization and magnetic polarization evaluated 
at the present. For a closed geometry, the integral is replaced by a sum over q/\K | = 3, 4, 5 . . . Note that there is no 
cross correlation Cf B or Cf B due to parity. 

We caution the reader that power spectra for the metric fluctuation sources Phil) — (h* (q)h(q)} must be defined in 
a similar fashion for consistency and choices between various authors differ by factors related to the curvature (see JlJ] 
for further discussion). To clarify this point, the initial power spectra of the metric fluctuations for a scale- invariant 
spectrum of scalar modes and minimal inflationary gravity wave modes B are 



P$ (q) cx 



1 



v(v 2 + 1) 



P h (?) a tanh(™/2) , (44) 

-*(v L + 1) 



where the normalization of the power spectrum comes from the underlying theory for the generation of the pertur- 
bations. This proportionality constant is related to the amplitude of the matter power spectrum on large scales or 
the energy density in long- wavelength gravitational waves |fl8|| . The vector perturbations have only decaying modes 
and so are only present in seeded models. The other initial conditions follow from detailed balance of the evolution 
equations and gauge transformations (see Appendix |a|) . 

Our conventions for the moments also differ from those in 0,^5) . They are related to those of |lj] byQ 



1 Footnote 3 of M incorrectly gives the relation between Q and A T . 
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FIG. 1. The scalar (left) and tensor (right) angular power spectra for anisotropies in a critical density model (thick lines) 
and an open model (thin lines) with flo — 0.4. Solid lines are Cf e , dashed Cf E and dotted Cf B . 



(2/ + l)A^=e[ 0) /(2.) 3 / 2 , 
(2i + l)Ag> = V2epV(27r) 3/2 . 



(45) 



where the factor of v2 in the latter comes from the quadrature sum over equal m — 2 and —2 contributions. 



Similar relations for A^'^L occur but with an extra minus sign so that Cc,i — —Cf E with the other power spectra 



*(E,B) 



unchanged. The output of CMBFAST continues to be Cc,i with the sign convention of |lif| . In the notation of jl5| 



the temperature power spectra agree but for polarization C. 



EE,BB 



Cf ' C /2 and Cf E = -Cj G l^. 



IV. RESULTS 



We now employ the formalism developed here to calculate the scalar and tensor temperature and polarization power 
spectra for two CDM models one with critical density and one with f^o = 1 — = 0.4 with initial conditions given 
by Eq. ([5]). In general, there are two classes of effects: the geometrical and dynamical aspects of curvature. 

On intermediate to small scales (large I), only geometrical aspects of curvature affect the spectra. Changes in the 
angular diameter distance to last scattering move features in the low-f2o models to smaller angular scales (higher I) 
as discussed in 



III 



Since the low-/ tail of the .E-mode polarization is growing rapidly with I, shifting the features 
to higher I results in smaller large-angle polarization in an open model for both scalar and tensor anisotropies. The 
suppression is larger in the case of scalars than tensors since the low-Z slope is steeper 

The presence of curvature also affects the late-time dynamics and initial power spectra. As is well known, the scalar 
temperature power spectrum exhibits an enhancement of power at low multipoles due to the integrated Sachs- Wolfe 
(ISW) effect during curvature domination. This does not affect the polarization, assuming no reionization, as it is 
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notation as Fig. |l] (thick: flat; thin open). Dotted lines represent negative correlation. 
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generated at last scattering. However it does affect the temperature-polarization cross correlation (see Fig. ||). In an 
open universe, the largest scales (lowest I) pick up unequal-time correlations with the ISW contributions which are of 
opposite sign to the ordinary Sachs- Wolfe contribution. This reverses the sign of the correlation and formally violates 
the predictions of Jl9| . In practice this effect is unobservable due to the smallness of signal. Even minimal amounts 
of reionization will destroy this effect. 

Open universe modifications to the initial power spectrum are potentially observable in the large angle CMB 
spectrum. Unfortunately subtle differences in the temperature power spectrum can be lost in cosmic variance. While 
polarization provides extra information, in the absence of late reionization the large-angle polarization is largely 
a projection of small-scale fluctuations. Nonetheless in our universe (where reionization occured before redshift 
z 1=3 5) the large-angle polarization is sensitive to the primordial power spectrum at the curvature scale. Thus if the 
fluctuations which gave rise to the large-scale structure and CMB anisotropy in our universe were generated by an 
open inflationary scenario based on bubble nucleation, a study of the large-angle polarization can in principle teach 
us about the initial nucleation event fisj] . 

In summary, we have completed the formalism for calculating and interpreting temperature and polarization 
anisotropics in linear theory from arbitrary metric fluctuations in an FRW universe. The results presented here 
are new for non-flat vector and tensor (polarization) perturbations and we have calculated the scalar and tensor 
temperature and polarization contributions for open inflationary spectra. The open tensor perturbation equations 
have been added to CMBFAST which is now publically available. 

Acknowledgments: We thank the Aspen Center for Physics where a portion of this work was completed. W.H was 
supported by the W.M. Keck Foundation and NSF PHY-9513835; M.Z. by NASA Grant NAG5-2816. 

APPENDIX A: EINSTEIN EQUATIONS 



In this Appendix, we complete the Boltzmann equations of § III by giving the Einstein equations for the metric 
and the matter. We begin with the background evolution and then proceed to the fluctuations. It is occasionally 
convenient to shift between different representations or gauges and thus we first discuss the transformations that link 
them. We then derive and present the Einstein equations for scalar, vector and tensor perturbations in a universe 
with constant comoving curvature in the synchronous and Newtonian gauges (see also (§)). 

1. Background Evolution 

The Einstein equations G^ v = SttGT^ express the metric evolution in terms of the matter sources. The background 
evolution equations are 

££ + 3(l + u; / )- = 0, 
Pf a 

4> + 2-<j) + a 2 V ^ = 0, (Al) 
a 

for the fluid and scalar field components respectively and 



2 



8nG 



n , +K =— a (Pf +P<P + Pv), ( A2) 



where Wf = Pf / Pf and P(j,{4>) was given in Eq. ( |l7| ) and p v = SHqCI/^/SttG is the vacuum energ 



2. Gauge Transformations 

To represent the perturbations we must make a gauge choice. A gauge transformation is a change in the correspon- 
dence between the perturbation and the background represented by the coordinate shifts 
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%i = x i ~\~ LQ 



(m) (A3) 



T corresponds to a choice in time slicing and L a choice of spatial coordinates. Since scalar and vector quantities 
cannot be formed from tensor modes (m = ±2), no gauge freedom remains there. Under the condition that metric 
distances be invariant, they transform the metric as Q 

_ j\{m) _ rp _ <^rp 

a 

£(m) = B (m) + l + kT ^ 

M m) =H { r m) --L--T, 
3 a 

= + kL . (A4) 

The stress-energy perturbations in different gauges are similarly related by the gauge transformations 

6f=6f + 3(l + w f )-T, 

Spf — Spf + 3c 2 fPf(l + Wf)-T , 
v { p ] = v ( } m) + L , 

7rJ m) = 4 m) . (A5) 

Note that the anisotropic stress is gauge-invariant. Seed perturbations are also gauge-invariant to lowest order, 
whereas a scalar field transforms as 

5~<f> = 5<t>-(f>T. (A6) 
The relation between the synchronous and Newtonian gauge equations follow from these relations. 



3. Scalar Einstein Equations 



With the form of the scalar metric and stress energy tensor given in Eqs. (A.4) and (|l5|), the "Poisson" equations 
become in the Newtonian gauge 



(fc 2 -3if)$ = AnGa 2 
fc 2 (* + $) = -8^Ga 2 ( W 4 0) +^ 0) ) 
and in the synchronous gauge 



(Pf6f + p.) + A(pf+Pf)v?>+vW]/k 



(fc 2 - 3K){h L + ^h T ) + 3^h L = AirGa 2 [p f 5 f + p s ] , 

h L + ^(l-3K/k 2 )h T = ~^Ga 2 [{pf+p f )vf ] +v^]/k : 



■■ a ■ 
h,L + -riL 
a 



a 



h T + -hr - k 2 (h L + -h T ) 



= -^Ga 2 [-p f 6f + Sp f + -p s + p s ] 
= -87rGa 2 [p / 7rf +7rf ]. 



(A7) 



(A8) 



Two out of four of the synchronous gauge equations are redundant. 

The corresponding evolution of the matter is given by covariant conservation of the stress energy tensor T^: 



Sf = —(1 + Wf)kVf — 3-Swj + Ss , 



(1 + Wf)vf } ] = -(1 + w f )-(l - 3w f )vf ] + w f k 
1 J a J 



6 Pf /p f -l(l-3K/k 2 )nf ] 



c(o) 



(A9) 
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for the fluid part (see e.g. [pOfl). The gravitational sources are 

S 6 = -3(1 + «>/)$, = (l + u>/)k* 



-3(l + u>/)/i£, Si 



(0) 



0, 



(Newtonian), 
(synchronous) 



(A10) 



These equations remain true for each fluid individually in the absence of momentum exchange, e.g. for the cold dark 
matter. The baryons have an additional term to the Euler equation due to momentum exchange from Compton 
scattering with the photons. For a given velocity perturbation the momentum density ratio between the two fluids is 

R _ Pb + Pb _ 3ps 

A comparison with photon Euler equation ([33]; 1 = 1) gives the source modification for the baryon Euler equation 



(AH) 



For a seed source, the conservation equations become 



(A12) 



6 (0) = _ 4 ^(o) + k 
a 



independent of gauge since the metric fluctuations produce higher order terms. 
Finally for a scalar field, tp = <j> + 8<f), the conservation equations become 



8<t> + 2 -86- 
a 



where 



are the gravitational sources. 



Set, = 




(fc 2 + a 2 V^)5(j) = S 



(Newtonian), 
(synchronous), 



(A13) 



(A14) 



(A15) 



4. Vector Einstein Equations 

The vector metric source evolution is similarly constructed from a "Poisson" equation: in the generalized Newtonian 
gauge 

V + 2-V = -SirG^ipfir^ + w^/k , (A16) 
a 1 

and for the synchronous gauge, 

h v + 2-h v = -STvG^lpf^P + n^)/k 2 . (A17) 
Likewise momentum conservation implies the Euler equation 

= -(1 - 3c})-vf } ~ - 2K/k^f + fiCt), (A18) 

1 1 a 1 2 1 + Wf 1 

where recall ci =Pf/pf is the sound speed and the gravitational sources are 



5(1) = J V+ (1 - 3c /)~ y ' (Newtonian), ^ A19 ^ 
[ , (synchronous) . 

The seed Euler equation is given by 

{,« = _4-«« - ifc(l - 2K/k 2 )^. (A20) 

Again, the first of these equations remains true for each fluid individually save for momentum exchange terms. The 
baryon Euler equation has an additional term in the source of the same form as Eq. (A.12) with m = — * m = 1. 
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5. Tensor Einstein Equations 



The Einstein equations tell us that the tensor metric source is governed by 



. a 



H + 2-H + (k 2 + 2K)H = %TiGa 2 [p f Tif ] + ir {2) ] , (A21) 



for all gauges. 



a 



APPENDIX B: RADIAL FUNCTIONS 



It is often useful to represent the eigenmodes in a spherical coordinate system (%, 9, (f>) , where x is the radial 
coordinate scaled to the curvature radius. Here we explicitly write down the forms and properties of the radial modes 
in an open geometry and describe the modifications necessary to treat closed geometries. 

By separation of variables in the Laplacian, we can write 

*G? = £H) VM2Z + 1) s a^ ( X , v) e YT(n), (Bl) 
i 

and the goal is to find explicit expressions for s a, . Here the ^-weights are set to reproduce the flat-space conventions 
of spherical Bessel functions (see also jlj). We proceed by analyzing the lowest j — min(|s|, \m\) harmonic 

—n . . . n '^ h ,, A]m] , 
±2 G r i L oc (mi ± i^tirht ± im 2 f 2 Q^l , (B2) 

where mi and m 2 form a right-handed orthonormal basis with n. We can now determine s a, from the radial 
representation of Q( m ) plj 



^ 13 (x,^) = w|^c S chx$r(x), 



^^'Vl ^X'+i) ^^' (B3) 



r (mm) i (mm) ■ -i i r (2m) (m) ■ - n(m) 

tor a v — <Pi ; similarly tor ±2 a\ = ± ^p ^ v , 



£ ' 1)(X ' V) = W (J + 4)(^+l) CSChX [COthX$r(x) + <f>r ' (x)] 



£ ' 2)(X ' V) = + ^ + [*r"(x) + 4coth X $r'(x) - (v 2 1 - 2coth 2 X ) $ r(x) ] 



(B4) 



and 



A (0) (x^) = o, 



^ 2)(X ' V) = ^+4) (,2 + 1) [<1>r(x) + 2c ° th * $r( * )] ' 



(B5) 
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for m > 0. For m < 0, /3; = — /3| m-1 while the other two functions remain the same. Here (\) is the hyperspherical 
Bessel function whose properties are discussed extensively by ||. 

The overall normalization of the modes here has been altered from those of J(|m in the case of vector and tensor 
temperature modes such that 



?( m ) 



(0,") 



1 



21 + 1 



(B6) 



where the difference lies in the lack of curvature dependence in the relation. Our choice simplifies the equations since 
it preserves the flat space form of the equations locally around the origin. It also defines the normalization of the 
polarization modes with respect to through Eq. (B2). 

The properties of the hyperspherical Bessel functions imply useful properties for the radial functions. For our 
purposes, the important relations they obey are: 

d 1 
d X 21 + 1 



lyM+P^ -(1 + 1)^2 + (Z + 1)2$^ 



cothx^r = 



1 



~ l 21 + 1 

which define the series in terms of its first member 



(B7) 



v sinh x 

Notice that limjf^o ®i(x) = 3l(kr). 

From the recursion relations of one establishes the corresponding relations for the radial function 



d X ls 1 ' 21 
for the lowest j, where recall 



s u l-l 



s^l + l 



„(i m ) 



1(1 + 1) 



Ajm) 



(B8) 



(B9) 



S™1 



(l 2 - m 2 )(P - s 2 ) 







I 2 ' 




n 









(BIO) 



The construction of the higher S G" 1 via the recursion relation of Eq. ( |25| ) also returns the higher radial harmonics. 
A few useful ones are 



1 



1 



*r(x), 



{u 2 +4) (1/2 + 1) 



(Bll) 



Furthermore, the recursion relation obeyed by the higher radial harmonics is the same as Eq. (B9), by virtue of 
Eq. ( p5| ) and explicit substitution of the radial form Eq. ( |C3| ) . This j-independence of the recursion relation implies 
that <pf m ^ is a sol ution to the temperature hierarchy Eq. ([33|) for any j and aids in the construction of the integral 
solutions in § QID . 

Finally, the radial functions for a closed geometry follow by replacing all v 2 + n, where n is integer, with v 2 — n 



and trigonometric functions with hyperbolic trigonometric functions (see |p|j21[| for details) 



APPENDIX C: DERIVATION OF THE NORMAL MODES 



We would like to describe the spatial and angular dependence of the normal modes s GJ n (a;, n) in a coordinate- free 
way by constructing them out of covariant derivatives of Q( m ) contracted with some orthonormal basis (n, rhx, 7712)- 
The lowest j = max(|m|, \s\) modes can be written as P|§], 
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±2 G? <x (mi ± irn 2 y (mi ± imn) h Q^ , (CI) 



and satisfy (Appendix O) , 



s GT(x,n) = (- i ) i y^I[ s y™(n)]exp^(f,fc)] 



(C2) 



with / = j. We demand that the higher /-modes also do so, to maintain the division of spin and orbital angular 
momentum defined in flat space pj]. 

We begin the construction by choosing some arbitrary point xq, and using a spherical coordinate system around it, 
x — xq — \J — K x(— n). Now n defines both the intrinsic angular coordinate system and the angular coordinates for 
the spatial location x(x,n). This reduction in the dimension of the space is sufficient since the end goal is to derive 
how the intrinsic and orbital angular dependence in the same direction h adds. In physical terms, only those photons 
directed toward the observer can contribute to the local angular dependence there. First expand the lowest mode in 
spin-spherical harmonics 

s Gf(x,n;u) = J2H) l V^l+T) Jr\x,y) ,17» , (C3) 
i 

where recall that the dimensionless wavenumber is v = q/y/—K. We obtain the explicit expressions for s ap m ^ and 



7 

their recursion relations in Appendix M by simple comparison between equations ( |Cl| ) and ( |C3| ) . At the origin they 
satisfy 



s a l 



(«) 



which both fixes the normalization of the modes and manifestly obeys Eq. flC2|) . As x —> 0, only the local angular 
dependence remains, as expressed in the Kronecker delta of Eq. (|C4| ). Because the spatial variation of the normal 
mode Q( m ) across a shell at fixed radius x must be added to the local dependence, even a mode of fixed j has a sum 
over all I in its angular dependence which contributes at any other point. 

This generation of higher /-structure as X increases suggests that we can use the radial structure of S G™ to generate 

the higher /-modes. From the radial recursion relation for gOif Eq. (B9), let us make the ansatz 



I V \JTTLS 

-.n l ( s G™)|i = [ s k™ s G™j - ,5/4+1 s^I+i] _ *777TTi s^T" 



V=I v * * y| * 21 + 1 13 1 s l - 1 s 1+1 s 1+iJ 1(1+1) 
That this series generates modes with the desired properties can be shown by returning to the spherical coordinate 



system. By explicit substitution of the radial form for S G™ of Eq. (C3) and by noting that in this coordinate system 



we obtain 



' AsGT)\i = -4risGT), (C6) 



•^'"•''^^'''^aTT 1 ^^ 1, (C7) 



(up to a phase factor) as desired. Since we have shown this for an arbitrary point, it is clear that Eq. (C2) holds in 
general. Note that this construction requires 

/ ^\LGT 1 i rLGT 2 2 ]\ = ^jS ll , h S„ 2 , (C8) 

for all x, as in the flat case of Eq. (E4), and defines our normalization convention. 
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